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Short Gamma Ray Bursts (SGRB) are believed to originate from the merger of two compact objects. If this 
scenario is correct, SGRB will be accompanied by the emission of strong gravitational waves, detectable by 
current or planned GW detectors, such as LIGO and Virgo. No detection of a gravitational wave has been made 
up to date. In this paper I will use a set of SGRB with observed redshifts to fit a model describing the cumulative 
number of SGRB as a function of redshift, to determine the rate of such merger events in the nearby universe. 
These estimations will be used to make probability statements about detecting a gravitational wave associated 
with a short gamma ray burst during the latest science run of LIGOA^irgo. Chance estimations for the enhanced 
and advanced detectors will also be made, and a comparison between the rates deduced from this work will be 
compared to the existing literature. 

I. INTRODUCTION 

Gamma Ray Bursts (ORB) are intensive bursts of high-energy gamma rays, distributed uniformly over the sky, lasting millisec- 
onds to hundreds of seconds. Several thousands of bursts has been discovered to date, with the very prominent feature of a 
bimodal distribution of the durations of the bursts, with a minimum around 2 seconds (B O. Bursts with a duration shorter 
than 2 seconds are called short GRB's, and bursts lasting longer than 2 seconds are labelled long GRB. Long GRB's have been 
associated with star-forming galaxies OIHO and are believed to be created by stellar core-collapse. In fact, several long GRB's 
have been associated with observed supernovae 161171 l8ll9HT0l. 

The origin of shorts GRB's, on the other hand, had been a mystery for a long time, although some time ago it has been proposed 
that they arise from the merger of two compact objects, like a binary neutron star (BNS) or a neutron star-black hole binary 
(BHNS)|[TT||T2]|. This picture for short GRB's is mostly accepted today, especially strengthened by recent observations of the 
SWIFT satellite |TT, 14,15]: 

First, while some of SGRB occur within starforming galaxies (like GRB050709), other SGRB's are being found in the outskirts 
of elliptical galaxies, without ongoing star formation (e.g. GRB050509B |[T6ll ). This can be explained by the assumption that 
the merger gets a high kick velocity when one of the components undergoes a supernova explosion. Then, due to the long time 
scale of 10-1000 Myr before the two components will merge, the system has enough time to leave the host galaxy and create 
a (S)GRB outside, which is expected to happen for a substantial fraction of mergers liTTl . Also, given the long time scale, the 
galaxies themselves have enough time to evolve to a late-type galaxy. 

Another support for the merger model comes from SWIFT observations of weak afterglows for some short GRB's |[T8l[T9ll20]| , 
which are much weaker than that of long GRB's. The weaker afterglow can be explained by the lower energy emitted in a 
SGRB, but also by an environment with a much lower density surrounding the GRB. This is consistent with both the merger 
hypothesis ||2T]|22| and the assumption that afterglow radiation is created by external shocks. This evidence favors the model that 
SGRB's originate from the merger of two massive objects, taking place in the outskirts of evolved galaxies, with the mechanism 
for creating the electromagnetically radiation the same as assumed for long GRB's (internal shock model) [23^, 24 1 . Therefore, 
the time- scale between any outgoing gravitational wave and the onset of the electromagnetic radiation should not exceed some 
milliseconds. 

Besides the merger scenario it is assumed that some of the short GRB's are caused by soft gamma repeaters (SGR), fast rotating 
magnetically-powered neutron stars, creating 'star quakes' in the crust from time to time and generating bursts of gamma radia- 
tion ||25l[26l. It has been estimated that up to --25 % of all SGRB's are caused by SGR's f27lE8|. If a reasonable large fraction 
of SGRB's is indeed created by the merger of two compact objects, they represent putative sources of gravitational waves and 
could produce a measurable signal in current or planned gravitational wave (GW) detectors, such as LIGO and Virgo. 
The LIGO detectors, described in detail in [29 , 30), consist of three kilometer-length, orthogonal interferometers at two sites. 
One detector is located at Hanford, WA (US) and the other at Livingston, LA (US). Both detectors contain a 4 km long inter- 
ferometer, while the Hanford detector additionally contain a 2 km interferometer housed in the same vacuum tube. The Virgo 
detector 1 3 1 1 is located near Pisa in Italy and consists of a 3 km long interferometer. All detectors work at or close to their design 
sensitivity and the LIGO detectors recently finished a 2-year data-taking run (S5) from November 2005 to October 2007. Virgo 
joined this data- taking effort in May 2007. Several results on searches for a merger signal have been published on data taken in 
earlier science runs |32, 33 , 34, 35 , 36] and analysis is finishing for the recent S5 run. An analysis using S5 data for a search of 
GW associated with GRB 070201 has been performed 1371 |38l, but so far, no GW has been detected. See 1391 for an overview 
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FIG. 1: Illustration of the different rate functions used for the models which are fitted to the observed data (see ||4T]|43l for details). Each 
model is scaled to unity at zero redshift. 



on gravitational wave astronomy and the status of the facilities. 

This paper will first describe the model fitted to the observed data, before the data and the fit results are presented. A comparison 
with rate estimations from pulsar observations and population synthesis follows, before chance estimations for the detection of 
GW's in current and planned detectors will be given. 



II. DESCRIPTION OF THE MODEL 



The aim of this section is to describe a most general model for the rate of astronomical objects, as a function of their cosmological 
redshifts. These models follow descriptions by Chapman 001 |4T1 and Guetta f42l and will be used to model the distribution 
of short GRB's in the following section. This section does not contain any new outcomes, but for clarity and for completeness 
every piece of this model is shown in details. 

This model, stating the number of objects with a redshift smaller than some redshift z*, is: 



Ni..)o.r,.^'^ HL)dL. (1) 

Jo 1 + ^ dz Jl ■ (P,- z) 

In this equation N{z^) is the number of SGRB with a redshift smaller than z^, R{z) is the rate-function (in units per volume) 
at a redshift z, ^{L) is the luminosity function and dV{z)/dz is the volume of a co-moving shell at redshift z. The integration 
over the luminosity function is taken out between a minimum luminosity Lmin{z), determined by the distance of the source and 



the satellite threshold, and an constant upper bound Lmax for computational reasons (see Section [irB] for details). 
The rate function R{z) describes the change of the intrinsic rate of objects as a function of redshift z using several approaches: a 
function that follows the star-formation rate, and two functions following a delayed star-formation rate. They will be described 
in Section irTAl in more detail. 

The luminosity function ^{L) describes the distribution of sources as a function of their luminosities; this could follow a single 



power-function, a Schechter function or a log-normal distribution. These functions are described in more detail in Section II B 
with a detailed derivation of the norming of these functions in Appendix |A| These functions have one to three free parameters, 
which are the ones being fitted. 



A. The rate function 



This section describes the rate functions which are used to fit the model to the data given in eq. ([T]). I will follow the same 
notation and enumeration as used in [ 42il : 
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1. A rate that follows the star formation rate (model "SF2" as described in 142114311 ): 

^/ N / N 23 6X0 (3.4 

R{z) = RSF2{Z) = Rs,o (2) 

exp (3.4 + 22 

2. A rate following the merger rate of two compact objects, as derived in |42| from six observed double neutron stars ll44ll . 
This rate is following a time-delay distribution (ex 1/r): 

R{z) = Rm,o / dr RsF2 {z{t - r)) /r (3) 
Jo 

3. A similar rate that follows SF2 with a constant time-delay distribution : 

R[z) = Rd,o / dr RsF2 {z(t - r)) . (4) 

^0 

In these equations t{z) is the look-back time for a redshift z and z{t) the inverse of that. Fig. [TJcompares these rate functions up 
to a redshift of 2 (scaled to be unity for zero redshift). 

B. The luminosity function 

The luminosity function describes the distribution of luminosities of the sources, for which often a power-law function or a 
schechter function is used in astrophysical context. These functions define the percentage r] of the sources detectable by detectors 
such as HETE-II and SWIFT: 

ri = ^ f^^^ HL) dL (5) 




FIG. 2: Plot of the approximate luminosity rate against the measured redshift for all the short GRB from Table [T| as blue circles. Also shown 
(solid red line) is the detector threshold, assuming a required minimum flux of 10~^ ergs /cm^ to produce a trigger (see Figure 12 in f45\). 
The area below the solid curve is not accessible by the instruments, and consistently all the GRB's lie above this line. The most left dot in this 
plot is GRB 080121, located at a very small redshift of 0.046 with a consistent luminosity 1461 . 
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This function describes the fraction of the sources which can be seen by a detector, given the luminosity function ^{L) = dN/ dL 
and a norming constant The lower limit of this integral depends on the detector threshold, which is ~ 1 photon/cm^ /sec for 
both HETE-II or SWIFT ||42l|45l, or roughly ~ l^~^ergs/ cm? (see Figure 12 in B31 ). The relationship for the minimum 
luminosity is simply 

47r 2 

Lniin{z) Anm(^) (6) 

Given the redshift of the source the corresponding minimum luminosity can be calculated, which acts as the lower integration 
limit in eqs. ^ and The upper bound for this integration, Lmax, is set to 10^^ ergs, which has been shown to be a reasonable 
value - using a larger upper limit has little effect on the outcome of the fits. 

The following list summarizes the different luminosity functions that are considered in the fitting model ([T]): 

1. A single power law distribution, which is often used to describe the pdf of luminosity in astrophysics (two parameters: 
and a): 

HL) = $0 (1^) (7) 

2. A broken power law distribution, describing e.g. two underlying populations in the luminosity [42J (four parameters: <l>o, 
I/O, a and P)\ 

^(L)= $0 (^)"" forL<Lo (8) 
^(L)= $0 (^)"^ forL>=Lo (9) 

3. The Schechter distribution, as used for example in ref |47 1 (three parameters: and a)\ 

HL) = $0 (^^) exp(-L/Lo) (10) 

4. A log-normal distribution, describing a standard candle, e.g. a population with about the same luminosity (following 1411 . 
with three parameters: ^o, Lq and a): 



The constant is always chosen so that the integral between the lower limit Lv and the upper limit for the luminosity is 
unity: J^^ ^{L)dL = 1. The used integration limits are Lv = 10^^ ergs and = L^ax = 10^^ ergs. These integration limits 
have been chosen because the derived luminosities for all used SGRB's given in table [l|range between 10^^ ergs and 10^^ ergs, 
well within the integration limits. It also has been shown that the results of the fit do not change significantly when widening the 
integration limits in either direction. The derivation of the constants $o can be found in Appendix [A[ 



III. USED DATA AND FIT RESULTS 



A. Used data set 



The data used to perform the fit consist of 11 short GRB's |75 | with well-defined redshift measurement up to z=0.5, detected 
by various satellite missions between May 2005 and September 2008. Some properties of these GRB's are given in Table |l] 
This table include fluences and luminosities for the GRB's as well, which are not used in the fitting procedure, because these 
values differ much between the instruments and the energy ranges; they are only used for a consistency check (see below). The 
cumulative number of the GRB's in Table[l|is used to fit the model given in eq. ([T]). 

Figure |2] shows the distribution of the luminosity rates (the luminosity divided by the Tgo time) of the SGRB's and the solid red 
line corresponds to the approximate threshold of the satellites (using a ballpark threshold of 10~^ ergs / cm'^ /s). As expected 
(and for a sanity check), all GRB's lie above the threshold line, including GRB 080121 (the most left dot in this plot) with 
a redshift of only 0.046. This particular short GRB has a very low luminosity, several orders of magnitude less than typical 
short-hard bursts [46 J . but is consistent with its distance. 
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FIG. 3: Histogram of the distribution of rates from the 24 accepted fit results, as described in the text. The rate values for this plot are also 
listed in the last column of Table HH 



B. Results of the fit 



The fit of the model was performed by a least-squares fitting function from the scipy module of python 11591 , which uses a 
modified version of the Levenberg-Marquardt algorithm to minimize a given function, which is the difference between the 
observed data and the fit (eq. ([T])). This fit depen ds on some parameters: one overall normalization factor and one to three 
parameters used in the luminosity function (see Sec. IIB ). Given these fit-parameters the -value of the fit is calculated, as well 
as the Kolmogorov-Smirnov (KS) probability, computed by taking the largest difference between the cumulative histogram of 
measured and estimated datapoints. From this difference a probability a is derived, at which the null hypothesis (i.e. samples 
are drawn from the same distribution) is not rejected (see e.g. 1601 for details). 

Besides performing the fit using only one luminosity function, fits were performed including two different luminosity functions 
in the case of more than one underlying population (see e.g. |41 1 for a similar approach). This might be the case, because of the 
contribution of SGR's to short GRB's, which do not originate from a merger. 

To take into account only reasonable models, fits with a KS probability less than 95% or with a value per d.o.f. of more than 
5 have been rejected. Because an estimated fraction of ~25% of all SGRB's are caused by SGR's LIT^^SJ. the major population 



GRB 


redshift 


duration [s] 


fluence [IQ-^^] 


luminosity [10^^ erg] 


050509Bf24l 


0.226 


0.040 


0.095 


0.012 


0507091 24 1 


0.1606 


0.22 


10.0 


0.175 


0507241 24 1 


0.257 


1.31 


9.9 


1.67 


0605051 481 


0.089 


4.0 


6.2 


0.16 


061006149 501 


0.4377 


0.42 


14.2 


8.27 


061210150 51| 


0.4095 


0.19 


11 


5.52 


0702091 52 1 


0.314 


0.10 


0.22 


0.059 


0704061531 


0.11 


0.70 


0.36 


0.009 


0707241 54 1 


0.457 


0.40 


0.80 


0.517 


0712271551 


0.383 


1.8 


2.2 


0.933 


080121|46| 


0.046 


0.7 


1.0 


0.004 



TABLE I: Set of 1 1 short gamma ray bursts with certain redshift measurements and Tgo duration equal or less than 4 seconds, and with a 
redshift smaller than 0.5, used for the fit. The data for these GRB's are taken from |56 57 1 and from diverse GCN circulars (see |58 |). The 
fluence and the luminosity of these GRB's are listed as well, but this information is not used for the fitting itself in this work. The luminosity 
has been calculated using the fluence and the luminosity distance to each GRB. 



6 



Model 


first model 


second model 


KS 


Chi 


rate 






a, a 




Lhreak 










Myr-^L-^ 


SF2 power 


— 


2.0 


— 





— 


— 


98.7 


0.15 


2530 


SF2 power/power 


— 


2.0 


— 





2.2 


— 


99.1 


0.18 


1710 


SF2 power/schechter 


— 


2.1 


— 


49 A 


-16.9 


— 


99.8 


0.16 


3480 


SF2 power/broken 


— 


2.1 


— 


48.0 


1.1 


0.8 


99.8 


0.19 


3490 


SF2 schechter/schechter 


47.6 


0.2 


— 


50.3 


-0.2 


— 


96.4 


0.12 


670 


SF2 schechter/lognormal 


48.7 


-8.2 


— 


47.2 


1.8 


— 


97.2 


0.10 


1070 


SF2 schechter/broken 


48.5 


-24.4 


— 


46.8 


-0.5 


2.4 


98.6 


0.16 


570 


SF2 lognormal 


46.8 


14.7 


— 


— 


— 


— 


98.6 


0.17 


2370 


SF2 lognormal/lognonnal 


47.2 


1.6 


— 


50.6 


0.4 


— 


96.8 


0.10 


890 


SF2 lognormal/broken 


47.2 


1.6 


— 


48.3 


1.1 


0.5 


96.8 


0.13 


890 


SF2 broken 


45.8 


-4.2 


2.0 


— 


— 


— 


98.0 


0.20 


1640 


SF2 broken/broken 


47.2 


0.7 


2.7 


48.0 


1.3 


0.9 


97.7 


0.23 


430 


Merger power 


— 


2.2 


— 


— 


— 


— 


99.1 


0.14 


3160 


Merger power/power 


— 


2.3 


— 


— 


2.0 


— 


99.7 


0.18 


3370 


Merger power/broken 


— 


2.1 


— 


47.5 


-1.0 


8.5 


99.7 


0.20 


1560 


Merger schechter/schechter 


47.6 


0.4 


— 


50.5 


0.2 


— 


96.2 


0.13 


590 


Merger schechter/broken 


49.4 


-5.0 


— 


46.9 


-0.9 


2.7 


96.6 


0.19 


240 


Merger lognormal/broken 


46.9 


-1.6 


— 


46.0 


-0.2 


1.4 


98.1 


0.14 


890 


Merger broken 


46.0 


-13.0 


2.3 


— 


— 


— 


97.7 


0.21 


1230 


Delay power 




2.3 










95.2 


0.18 


2740 


Delay power/power 




2.6 






2.3 




96.8 


0.22 


1770 


Delay power/broken 




2.2 




47.4 


-0.4 


8.6 


97.6 


0.27 


1450 


Delay schechter/schechter 


47.6 


0.5 




51.0 


0.2 




96.7 


0.14 


590 


Delay lognormal/broken 


46.8 


1.7 




48.0 


1.3 


0.6 


95.5 


0.12 


830 



TABLE II: Table containing the 24 selected fit results with a KS -probability of at least 95% and a P^r d.o.f. less than 5. These 24 models 
are being used for the further analysis. The models in bold font have the largest and smallest rates, and the value for Lq are given in units of 



of SGRB's is likely to come from the merger of two compact objects. Therefore the model with the larger rate prediction is 
chosen for any two-luminosity-model fit. Table [Il| summarizes the final sample of models, including their fitted parameters and 
the goodness-of-fit values. It is not surprising that many models with the SF2 rate-model fit the data well, because it has been 
shown that the distribution of SGRB's should follow the SF2 rate pretty closely ifTTllMl . 

To obtain a rate estimation for a model, each model (i.e. equation ([T])) is integrated up to a distance of 100 Mpc (corresponding 
to a redshift of 0.02258), and rescaled by the number 11 of GRB's used in the fit. This number corresponds to the probability 
that a given GRB lies within 100 Mpc. This number must be further multiplied by a factor 15/28 800 x 100. 15 is the approximate 
number of short GRB's occurring during one year, 28 800 is the approximate number of Lio within 100 Mpc |76|, and f^^ = 
100 is the approximate beaming factor of a SGRB 1241 162|, to take into account all arbitrarily oriented mergers as well. The 
values of the rates calculated that way are given in the last column of Table [ll| 

Figure |3] shows the distribution of the rates of these 24 models, spreading over a large range of values. This is not a surprise 
when considering that the fit itself is performed up to a redshift of 0.5, but the rate is estimated within a redshift range of only 
z 0.03. It also should be noted that the outcome of the fits is very sensitive on the given input data; using only 10 out of the 
1 1 GRB's to fit a model can change the rate estimation by an order of magnitude. Plots of all accepted models and the data are 
shown in Figs. [5][8]in the Appendix. 

As a result of the analysis 20% of the models give a rate larger than 2750 Myr~^ L^q^, 50% of the models yield a rate of 
> 1450 Myr~^ L^q^, while 80% of all models predict a rate of > 670 Myr~^ L^q^. Since the 50% value corresponds to a median 
value, this value will be used in the next section to compare with rate estimations. 
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FIG. 4: Plot showing the detection probabihty as a function of detector range (left plot for the BNS scenario, right plot for the BHNS). The 
diversity in the detection probability in the y-direction comes from the different fits. The red curve signifies the detection probability for the 
'worst' and 'best' model, while the yellow area shows the range of probability predicted by 20-80% of the models. The vertical lines denote 
the probability areas as predicted by the models. See also Table [V] 



IV. COMPARISON WITH OTHER RATE ESTIMATES 



In this section a rough comparison is made between the rates deduced in this paper and rates estimated elsewhere. Two cases 
are being distinguished: The rate of mergers of two neutron stars (BNS, Binary Neutron Star) and the rate of mergers of a Black 
Hole and a Neutron Star (BHNS). For BNS the rate is deduced from known binary pulsars in our Milky Way and are expected 
to be realistically at 50 Myr"^ L^q\ although they could be as high as 500 Myr"^ Lj~o^ The rates predicted for BHNS 

are much more uncertain, and have been modeled using population synthesis. Realistic rates for BHNS lie at 2.0 Myr~^ L^q^, 
although they could be as high as 60 Myr~^ L^q^ 1641 . 

The models used in this work cover a rate range of 240 Myr~^L^Q^ to 3500 Myr~^L^Q^, with a median value at 
~ 1450 Myr~^ L^q^. These values are consistently larger than the previous rate estimates. 

A similar work on using GRB's to fit the local observed population has been done in |[62ll , quoting two rates of 130 Gpc~^yr~^ 
and 400 Gpc~^yr~^ , respectively, for two different models on the origin of binaries. This corresponds to a rate of 460 Myr~^ L^q^ 
and 1400 Myr~^ L^q^ ; the latter value surprisingly close to the rate of SGRB's deduced in this work. This might be a hint that the 
rates estimations in |[64l|65l[66l are underestimated, but given the large uncertainties involved in the fit procedures in this paper, 
more statistics is needed to investigate this hypothesis. Table III summarizes the median/realistic rates deduced from several 
methods, papers and this work. 





realistic/median rate 
[Myr-^L-^] 


this work 


1450 


Ref. 1 62 1 model I 


460 


Ref. 1 62 1 model II 


1400 


BNS rate | 67 1 


50 


BHNS rate | 64 1 


2.0 



TABLE III: Comparison of realistic rate estimations for BNS and BHNS with the median values of the rates deduced in this work, compared 
to the outcome of a similar investigation by Guetta |62|. The first three numbers are derived from fits to observed short GRB's, while the last 
two rates are deduced from binary pulsar observations and population synthesis. It is interesting to note that the rates estimated from observed 
redshifts (rows 1 to 3) are consistently larger than the ones estimated otherwise. 
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GRB 


detectors 


range [Mpc] 


p [%] (BNS) 


p [%] (BHNS) 


GRB051114 


H1H2 


20.1 


0.02 - 0.08 


0.12-0.47 


GRB051210 


H1H2 


21.8 


0.02-0.10 


0.15-0.60 


GRB051211 


HlLl 


33.0 


0.10-0.41 


0.51-2.07 


GRB060121 


HlLl 


5.3 


0.00 - 0.00 


0.00 - 0.01 


GRB060313 


H1H2 


21.1 


0.02 - 0.09 


0.13-0.54 


GRB060427B 


HlLl 


56.2 


0.49-1.98 


2.49 - 9.77 


GRB060429 


H1H2 


33.0 


0.08 - 0.34 


0.51-2.06 


GRB061006 


H1H2 


21.9 


0.02-0.10 


0.15-0.60 


GRB061201 


H1H2 


30.5 


0.07 - 0.27 


0.40-1.64 


GRB070201 


H1H2 


15.5 


0.01 - 0.03 


0.05 - 0.21 


GRB070707 


H1H2 


31.3 


0.07 - 0.29 


0.43 - 1.76 


GRB070714 


HlLl 


17.2 


0.01 - 0.06 


0.07 - 0.29 


GRB070729 


HlLl 


52.8 


0.41 - 1.64 


2.07-8.18 


GRB070809 


H1H2 


10.6 


0.00 - 0.01 


0.02 - 0.07 


GRB070923 


HlLl 


19.7 


0.02 - 0.09 


0.11 -0.44 



TABLE IV: 15 short GRB's detected during the fifth LIGO Science run and VSRl of Virgo, that have no reasonable redshift estimation, and 
for which data is available from at least two detectors (indicated in the second column). The third column shows the approximate detection 
range of the second most sensitive detector at this time, using the SNR thresholds as described and the known position of the source. The 
fourth/fifth column gives the probability range for each GRB for the BNS/BHNS case. The total probability of a detection during S5 lies in 
the 20/80 range of '-^1.9 to ^7.7% with a median value around ^4.1% for the BNS case. When using the 'best' model, the total probability 
increases to 9.7% (see Table|V]for details ond for the values of the BHNS case). 

V. CHANCE ESTIMATION FOR LIGO AND VIRGO 

This section discusses the probability of a detection of gravitational waves associated with a SORB during 85. To calculate this 
probability the reachable distance of a detector to a source must be calculated, given the known horizon distance of a detector at 
a given time. 

The horizon distance of a detector is the distance at which a merging system would create a signal in this detector with a 
signal-to-noise ratio of 8.0, if the system is optimally located and optimally oriented. For the approximate calculation of the 
probabilities of discovering a SORB within the analysis of 85 data (and for simplification) two cases are taken into account for 
the merger: a low-mass pair consisting of two 1.4 Mq neutron stars (BNS) and a high-mass pair consisting of a 1.4 Mq neutron 
star and a 10 Mq black hole (BHNS). Looking at the ranges of the inspiral analysis in S5 (see, e.g. | 68 1), the horizon distance 
for detecting the BNS system is about 30 Mpc for Hl/Ll and 15 Mpc for H2/V1, while for the BHNS case it is 50 Mpc and 
25 Mpc, respectively. In the following it is assumed that Virgo behaves similar to H2 with respect to the sensitivity. 
The actual search for gravitational waves of mergers from a short GRB is split up into two cases, depending on the available 
data, and utilizing two different thresholds. If data is available from the two most sensitive detectors HI and LI (L-case), the 
threshold is set for them to be 4.25, and in the case that one of the detectors is H2 or VI (H-case), the threshold for the least 
sensitive detector is set to 3.75. It is very unlikely that a real gravitational wave, producing a signal just above these thresholds, 
will be identified as such, since it will be buried well within background noise. Therefore I will assume a rather conservative 
threshold for signal detection 50% above the actual used threshold. 

Taking this into account the distances to which a binary system can be detected for the L-case are ~35 Mpc (BNS) and ~65 
Mpc (BHNS) and for the H-case they are ~20 Mpc (BNS) and ~35 Mpc (BHNS), respectively. These are, however, only rough 
estimates of the distance which can vary up to ~ 20% because of fluctuations in the taken data (see, e.g. the Figure showing the 
reach of the detectors as a function of mass in | 68 1). 

Because a merging system is in general neither optimal located nor optimal oriented this has to be taken into account in form 
of the so-called antenna factors 1 69], which depends of the sky location of the putative source and the inclination angle l of the 
binary system. 

This inclination angle is the angle between the axis of total momentum of the system and the direction of the line of sight from 
the system to earth. For a system seen from earth face-on this angle is zero, while for a system seen edge-on this angle is 90°. 
Since the axis of total momentum of the system is the only significant direction and we are able to see the electromagnetic burst, 
taking place in a jet within some opening angle, it is a reasonable assumption that the inclination angle is smaller than or equal 
to the opening angle. For short GRB's, the opening angle was determined to be at about 0.16 radians 124 J . For the calculation of 
the probabilities it is assumed that the inclination angle is below ^ = 0.2. 
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Using the sky location it is possible to calculate the reachable distance for each GRB with the appropriate SNR threshold. The 
results of this calculation for all GRB's during S5 is shown in Table IV This table also contains the probabilities that a particular 
GRB lie within the detectable range as predicted by 20% (80%) of the models. The total detection probability that a gravitational 
wave will be seen in S5 data is 7.7% for 20% of the models, while 80% of the models predict a probability of at least 1.9%. The 
median probability lies at 4.1% (see also TablejV]). The probabilities for the BHNS case lie in the 20/80 range of 9.9% to 34.7%, 
with a maximum value of 41.9%, and a median probability of 20.2%. 



VI. OUTLOOK FOR ADVANCED DETECTORS 



Currently LIGO and Virgo are in the process of upgrading, and a science run with the enhanced detectors will start in summer 
2009. These detectors will be twice as sensitive than initial detectors fTOll and hence can observe an eight fold larger volume 
in the universe. For advanced LIGO, supposed to start data- taking in 2014, the increase in sensitivity with respect to the initial 
detectors is a factor of 10 iTOll , corresponding to an 1000-fold increase in rate. 

The chances are very good that during the runs of the enhanced and advanced detectors enough satellites are operating to deliver 
gamma-ray triggers. Currently three dedicated GRB missions are orbiting the earth, HETE-II|71|, SWIFTI 14. .151 and Fermi 
(formerly GLAST) 1721 . While HETE-II might last some more years, the SWIFT satellite is expected to operate for the next five 
years, covering at least the science runs with the enhanced detectors; SWIFT may be even operational into the runs with advanced 
detectors. Also the recently started Fermi satellite, launched in June 2008, is expected to work at least 10 years |73|, enough 
time to cover data taking periods with advanced LIGO and advanced Virgo. Last to mention is the Interplanetary Network (IPN) 
11741 , a network of spacecrafts equipped with gamma-ray detectors, such as Wind, NEAR, KONUS and Mars Odyssey. By 
measuring the time-delays of a trigger between these spacecraft it is possible to triangulate the approximate position of a GRB. 
These positions are not very precise in contrast to dedicated GRB's missions, but in most cases precise enough for an associated 
search of GW from this event. 

To estimate the probability of a detection with these detectors, a 2-year data taking run is simulated for each the enhanced and the 
advanced LIGOA^irgo detector network. 30 short GRB's have been simulated during this time, and the chance that at least two 
detectors will take data at the time of each GRB, assuming a duty-cycle of 80% for a single detector. The randomly selected sky 
location defines the sensitivity to all available detectors, of which the second most sensitive is being chosen, since this detector 
will limit the reach of the observation. Then for each model (taken from table [II]) the chance is being calculated that at least 
one of the 30 SGRB's is within detectable range and at least two detectors are 'data taking'. If that is the case, a 'detection' has 
been made. This procedure is repeated 5000 times for each model to estimate the chance to have a detection with a given model. 
This large number of trials ensures a relatively small Monte Carlo error of about 5%; the uncertainties in the duty-cycle and the 
number of SGRB per year have a much larger effect on the uncertainties. A change in the duty-cycle between 70% to 90% or a 
change in the number of SGRB's within 2 years between 25-35 yields an error of the MC results on the order of ~ 30% each. 
Figure]?] shows the outcome of this simulation for arbitrary sensitivities or detector ranges; the sensitivities for initial, enhanced 
and advanced detectors are marked by vertical lines. The spread in the detection probabilities come from the spread in the rate 
predictions of the different models and the yellow area shows the probability range predicted by 20-80% of the models. 
For initial LIGO (dotted vertical blue line), the predicted detection probabilities for the BNS case range from 0% (worst model) 
to 15% (most optimistic model). The 20/80-area covers a range between 3% and 13%, i.e. for 20% of the models the chance of 
a detection is at least ~13% and for 80% this chance is at least 3% (see also the second row in Table |V]). These predictions are 
consistently larger than the estimates in the first row in Table[V]using the real observed GRB's, because for most of S5 the Virgo 
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initial detectors (real GRB) 


0.7% 


1.9% 


4.1% 


7.7% 


9.7% 


3.7% 


9.9% 


20.2% 


34.7% 


41.9% 


initial detectors (sim GRB) 


1% 


3% 


7% 


13% 


15% 


7% 


17% 


30% 


51% 


59.6% 


enhanced detectors 


9% 


22% 


43% 


65% 


74% 


39% 


75% 


94% 


99.3% 


99.8% 


advanced detectors 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


100% 


100% 



TABLE V: The outcome of a Monte Carlo simulation (except the first line) to estimate the probability to have at least one short GRB within 
the reachable distance of the second most sensitive detector, given 30 SGRB's within the observing time of the network. A duty cycle of 80% 
for each single detector was used in the simulation. The first line was calculated using the properties of the known GRB's during S5. The 
numbers stated are the detection probabilities for the worst/best model, and the probabilities predicted by 20/50/80% of the models. For initial 
LIGO the predictions are not very promising, although the best model predicts a ^ 60% chance to see a detection for the BHNS case. The 
100 % values given for the advanced detectors are because of rounding effects; when doing sufficient simulations and considering the large 
uncertainties in the models they never reach exactly 100 % - although they are pretty close to it. 
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detector was not online. 

With the enhanced detectors the changes of detecting a gravitational wave associated with a GRB are much larger, and range 
between 22% to 65% (20/80-area) for the BNS case and up to 100% for the BHNS case. And for the advanced detectors finally, 
the Monte Carlo simulations predict a 100% chance for the detection of a gravitational wave, even for the worst model! The real 
chances for a detection are not exactly 100% (which is due to finite simulations), but very close. Then, at last, a new window to 
the universe should open up. Table|V]gives a summary of the detection estimations for initial, enhanced and advanced detectors. 
It has to be pointed out, as mentioned earlier, that these probabilities are based on a small data sample and therefore inherit large 
uncertainties. The largest uncertainties come from the small number input data itself, to which the fits are very sensitive to. Other 
uncertainties come from the assumed duty cycle of the detectors, the number of SGRB's detected by gamma-ray satellites during 
the science runs, and the SNR threshold of a putative signal. However, taking these uncertainties into account, the probability to 
have a detection of a gravitational wave associated with a short GRB with the enhanced detectors is about 50%, which is a very 
exciting prospects for the next science run expected to start in 2009. 

VII. SUMMARY AND CONCLUSION 

In this paper I made an attempt to estimate the local rate of short gamma ray bursts from a set of 1 1 short gamma ray bursts 
with reliable redshift measurements. I then compared them with values found in literature and used them to predict the chances 
of a detection of a gravitational wave assuming such events are mostly created by the merger of two compact objects. This 
assumption has gathered much support during the last years by observations with satellites like Swift. 

A model, describing the number of short GRB's as a function of redshift using different rate- and luminosity functions, has been 
used to fit the observed data. These fits, made with a modified version of the Levenberg-Marquardt algorithm, tried to minimize 
the difference between model and data, by altering the free parameters in the model. A KS test and a criterion was used to 
select only models with reliable good fits. Because these fits were made to a redshift of 0.5, but conclusions about the local rate 
are drawn from redshifts smaller than ~ 0.02, they have to be taken with large caution. Also, the variation of the input data has 
a large effect on the results and can change the results by an order of magnitude. 

Taking the median value of the fit results, which are spread over an order of magnitude, they look to be consistently larger than 
the predictions on the local rate made from binary pulsar observations and population synthesis. This is also the outcome in a 
similar investigation, leading to the possibility that the real local rate of mergers is larger than previously thought; but it is too 
early to draw firm conclusions on that. 

Although the rate estimations are relatively large, the prospects to discover a gravitational wave associated with a SGRB in the 
fifth science run in LIGO/Virgo are not promising; the probability for such a detection are well below 50%. Unlike larger are the 
chances to detect a gravitational wave in enhanced LIGOA^irgo: Monte Carlo simulations predicts a probability of ~ 50% for the 
BNS case and even ^90% for the BHNS case. For advanced detectors the chance of detecting a gravitational wave associated 
with a GRB turns out to be ~ 100%, even for disadvantaged models (for both BNS and BHNS cases). This large probability has 
to be taken with caution again, but even taking into account all the uncertainties they are very large. Such a coincident detection 
will not only give strong significance of the true nature of a real signal, it will also allow a much better distance estimation to 
cosmological sources, if the redshift of the GRB is known. Further observations of GRB associated gravitational waves therefore 
will not only increase the knowledge on stellar formation and galaxy evolution, but will have also a huge impact on the theories 
on large structures in the universe and on cosmology. 
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APPENDIX A: NORMING THE LUMINOSITY FUNCTIONS 

This appendix shows the derivation of the proper norming for the luminosity models. In the following the substitution x = L/Lq 
is used to calculate the norm 

1. A single power law: 

*o = (A1) 

With this norming it is obvious that the value of the integral J^^ ^{L)dL actually does not depend on the value Lq, but 
only on the exponent a. 



1 - 



1 



(A2) 



2. The broken power law: 

Integration of ^ gives the following normalization constant: 



' ^a+l _ a+1 /3+1 _ /3+1 ^ 

Xf) X\/ Xa X. 



-1 



and so the complete expression of the integration of the luminosity function from a lower bound xtox^ gives: 



if the value x lies within the interval [xv; xq] = [x\y; 1]. If x is larger than 1 the single power law can be used, scaled by a 
value K, so that the value at the intersection point (which is x = 1) falls together. Thus we require: 

^ox\^^^]=^^, ^ (A5) 



from which follows: 



/3+1 _ 6+1 

».*„£^_£^^ (AO) 

3. The Schechter and the log-normal distribution 

For the Schechter and the log-normal distribution the norming (i.e. the cumulative distribution) is computed numerically 
by numerical integration. 
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FIG. 5: Part I of the fits with the SF2 rate function. The x-axis shows the redshift and the y-axis the cumulative number of SGRB detected 
up to the given redshift. The red dots represent the observed data while the green line shows the best-fit model. In case of a two-model fit the 
thinner red and blue lines represent the contributions from both models. 
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FIG. 6: Part II of the fits with the SF2 rate function. The x-axis shows the redshift and the y-axis the cumulative number of SGRB detected 
up to the given redshift. The red dots represent the observed data while the green line shows the best-fit model. In case of a two-model fit the 
thinner red and blue lines represent the contributions from both models. 
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FIG. 7: Fits made with the delay rate function. The x-axis shows the redshift and the y-axis the cumulative number of SGRB detected up to 
the given redshift. The red dots represent the observed data while the green line shows the best-fit model. In case of a two-model fit the thinner 
red and blue lines represent the contributions from both models. 
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Rate-model: merger , Luminosity models: schecliter/schechter 




redshift 



FIG. 8: Fits made with the merger rate function. The x-axis shows the redshift and the y-axis the cumulative number of SGRB detected up to 
the given redshift. The red dots represent the observed data while the green line shows the best-fit model. In case of a two-model fit the thinner 
red and blue lines represent the contributions from both models. 
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